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Abstract. A one-dimensional chain of sporadic maps with asymmetric nearest 
neighbour couplings is numerically studied. It is shown that in the region of 
strong asymmetry the system becomes spatially fully synchronized, even in the 
thermodinamic limit, while the Lyapunov exponent is zero. For weak asymmetry the 
synchronization is no more complete, and the Lyapunov exponent becomes positive. 
In addition one has a clear relation between temporal and spatial chaos, i.e.: a positive 
effective Lyapunov exponent corresponds to a lack of synchronization and vice versa 
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An interesting property of extended systems is that they can exhibit complex 
behaviour both in space and time, that is the chaotic evolution of a spatial pattern. 
This kind of phenomena, called "spatio-temporal" chaos, has been the focus of 
considerable interest in different fields, as chemical react ion- diffusion systems, Benard 
convection, turbulence, modelling of brain fuction, and so on [1-2]. Recent experiments 
have revealed the emergence of large-scale spatio-temporal patterns of activity in 
many brain areas, such as the olfactory system or the visual cortex [3]. Coherent 
spatial patterns play also an important role in the behaviour of turbulent fluids. Here 
numerical simulations and experiments show the emergence of strong choerent vortex 
structures which are responsible for intermittency and possible anomalous dimensions 
in the scaling law [4]. Moreover choerent structures may play a non trivial role in 
the predictabilty problem based on two-dimensional turbulence, an important subjec 
from both theoretical and metereological point of view [5]. 

A direct study of the spatio-temporal behaviour from the full equations, such as 
the Navier-Stokes equations, in general, is quite complicated, even numerically. To 
overcome this problem, and gain more insight in the basic mechanisms, much work 
has been devoted to the study of models, such as coupled lattice maps. These are 
a crude, but nontrivial, approximation of extended systems with discrete space and 
time, but continuous states [2]. The simplest form is given by a set of N continuous 
variables Xi which evolve in (discrete) time as 

Xi(n + 1) = (1 - en - Pi) f{xi{n)) 

(1) 

+ aj(xi_i(n)) + (3if(x l+1 (n)). 

The index i denotes the i th site on a one-dimensional lattice, and only nearest 
neighbour interactions are taken into account. The strenght of the couplings and their 
spatial homogeneity, as well as the type of boundary conditions used, may depend on 
the the type of physical problem. For example in the case of shear flow, boundary 
layers or convection, there is a privileged direction. This can be introduced in the 
model (1) by taking asymmetric couplings [6,7]. 

Moreover in such systems the most appropriate boundary conditions are open 
boundary conditions. Thermal reservoir can be included by giving a special evolution 
law to the sites on the boundary. We do not consider this generalization here. 

In this letter we shall study a one-dimensional chain of coupled maps (1) with 
homogeneous asymmetric couplings and open boundary conditions. The equations of 
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motion are given by (1), with couplings 

a* =71, A =72 for % = 2, . . . , N - 1 

(2) 

ai =0, /3i=7 2 , «Ar = 7i, /3n = 0. 

Without loss of generality we assume 71 > 72- The system (1) with these couplings 
was studied in Refs. [7,8] for the case of chaotic single map f(x). In these papers 
the interest was for the possible emergence of spatial coherent patterns due to the 
asymmetrical couplings. 

A simple inspection reveals that the system (l)-(2) posses the uniform solution 

Xi(n) = x{n), x(n + 1) = f(x(n)). (3) 

To study the stability of the uniform state one can linearize eq. (1) about the uniform 
solution (3) and study the spectrum of the fluctuations [7] . This consists of a uniform 
eigenmode with eigenvalue po and N — 1 nonuniform eigenmodes with spectrum 

p(k) = po [1 - 71 - 72 + 2^7172 cos(fc)] (4) 

where k = irm/N with m = 1, 2, . . . , N — 1, and po = exp(Ao) where Ao is the largest 
Lyapunov exponent of the single map f(x). 

For 71 7^ 72 this spectrum posses a gap at k = since p(k — > 0) is less that the 
k = eigenvalue po- Therefore if 71 and 72 are such that po > 1 and \p(k)\ < 1 then all 
nonuniform fluctuations are stable. The only instability that is left is the instability 
to uniform fluctuations, which is inherent to the chaotic nature of the single map. 

^From this one concludes that the uniform state (3) is stable. This predictions for 
the case of a single chaotic map f(x) are essentially confirmed by numerical simulations 
[7,8]. The finite coherence length l c - that is: Xi(n) ~ x(n) for i < Z c , while for 
i > l c the xi are spatially irregular - is due to numerical noise. Indeed l c decreases 
logarithmically with the noise level in the numerical simulations. 

The open boundary condition can be seen as a defect in the chain. This scenario 
is essentially unchanged if one uses periodic boundary conditions and "softer" defects, 
such as, e.g., the interchange of 71 and 72 in a finite fraction of sites [8]. 

The above argument neglects the fluctuations of the chaotic degree along the 
trajectory, since the Lyapunov exponent Ao gives only the typical chaotic degree. If 
one considers a temporal window [t — At/2, t + At/ 2] it is possible to repeat the 
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arguments of Ref. [7] simply replacing Ao in equation (4) with the effective Lyapunov 
exponent [9] which measures the local exponential rate of growth for the tangent vector 
z around the time t: 

1 , |z(t + At/2)j 

We expect that the scenario discussed in Ref. [7] may fail if the fluctuations of XAt 
are strong. It is hence interesting to inquire what happens in the extreme cases. 

A rather natural way to address this problem is to study the system (l)-(2) with 
the single map 

f(x) = x + cx z modi (z>l). (6) 

For z > 2 the dynamical system x(n + 1) = f(x(n)) shows a sporadic behaviour 
[10], i.e. an initial disturbance grows in time as a stretched exponential with exponent 
less that one. The Lypunov exponent Ao hence vanishes. 

In the case of 1 < z < 2 the map (6) behaves as an ordinary chaotic system 
with positive Lyapunov exponent Ao- Here we obtain for the system (1) the scenario 
discussed in Refs. [7,8] for the logistic map. We thus find that for symmetric couplings 
the chain does not synchronize, while in the asymmetric case a finite coherence length 
may appear depending on the value of 71 — 72 and Ao- To be more specific we can 
distinguish two cases. In the first case po > 1 and \p(k)\ < 1 e.g. when z = 1.5, for 
which Ao = 0.56, 71 = 0.7 and 72 = 0.01 - so that all the non-uniform eigenmodes 
are stable an the system synchronizes with a finite l c [7,8]. In the second case - e.g. 
when z = 1.5, 71 = 0.7 and 72 = 0.1 - there exist some k for which \p(k)\ > 1, and 
the system does not synchronize. 

We now turn to the case of z > 2, when the single map is sporadic. The results 
we shall report are obtained from numerical simulations with 71 = 0.7 and 72 which 
varies in the range (0,71]. The typical sistem size considered is N = 200, but we 
checked some results also for larger values of N. 

If the couplings are symmetric the systems is not spatially synchronized, l c ~ 1. 
Moreover, even if the Lyapunov exponent of the single map Ao is zero, the Lyapunov 
exponent A of the global system is always positive confirming the non relevance of the 
uniform mode for the dynamics. 

As for the case of a chaotic map f(x), when the couplings are not equal two 
qualitative different behaviours appear depending on the value of 71 — 72. However, 



unlike the previous case, here the synchronization, when present, is complete, i.e. 
l c = N. In Fig.l we show r - the fraction of the time the system is completely 
synchronized - and A - the Lyapunov exponent - as functions of 72. 

A first regime, with r ~ 1, appears for 71 — 72 large, which with our parameters 
means 72 < 0.2. Here po = 1 and \p(k)\ is well smaller than one for all k, e.g., 
maxfc \p(k)\ = 0.85 for 72 = 0.2. According to the above argument the uniform solution 
is locally stable. Indeed we find that the system synchronizes with l c = N. A detailed 
analysis of the evolution reveals that the system spends most of the time completely 
synchronized. The ordered state is interrupted by intervals where the coherence length 
becomes very small. This is a very different qualitative behaviour with respect to the 
case z < 2 and the case discussed in Refs. [7,8]. In fact in these latter cases the 
system is always partially synchronized, with l c undergoing small fluctuations about 
its time average. All the simulations in this range of parameters gives a value of r 
very close to one. For example, when 72 = 0.1 we find, on a run of 5 x 10 6 steps, 
r = 0.997, while for 72 = 0.2, and on a run of 4.7 x 10 8 steps, we find r = 0.999. 
The Lyapunov exponent is always very small. It is difficult to decide if A is really 
going to zero or stays very small. This difficulty is due to the fact that every time the 
system desynchronizes, the local Lyapunov exponent increases. As a consequence the 
numerical value of A does not saturate in time but has kicks which increase its value, 
see Fig. 2. We note that this incertitude on the value of A is directly connected to that 
on t. 

A different behaviour appears for 71 — 72 small which, with our parameters, means 
72 > 0.35. In this regime po = 1 but \p(k)\ is larger than one for some values of k, 
e.g., when 72 = 0.35 we have max^ \p(k)\ = 1.04. Thus the uniform state is locally 
unstable. Indeed we find that the coherence length is almost always l c ~ 1 and the 
system is never completely synchronized - at least on a run of 3 x 10 7 steps. Moreover 
r is zero and A is positive. For example, A = 0.74 when 72 = 0.35. 

These two behaviours appear to be related in a continuous way through the 
crossover region where maxfc \p(k)\ is very close to one, for example max^ \p(k)\ = 0.92 
for 72 = 0.3. In this region A begins to be definitely non zero and r becomes smaller 
than one. For example, A ~ 0.05 and r ~ 0.9 again for 72 = 0.3 on a run of 3 x 10 8 
time steps. 

In order to stress the correlation between desynchronization and temporal 
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chaoticity, we show in Fig. 3 the effective Lyapunov exponent XAt(t) and l c (t) as 
functions of t. It is evident that a large XAtif) corresponds to Z c << N and viceversa. 
This happens both for small and large values of 71 — 72- 

We stress that since the map f(x) is sporadic, and hence does not posses an 
asymptotic measure, the initial condition is relevant. In fact while the qualitative 
behaviour described above does non change, at least for the values of 72 we used, the 
numerical value of A and r my depend on the chosen initial condition. 

This dependence is more dramatic in the crossover region, and makes difficult to 
define clear boundaries. 

We conclude with some remarks. The behaviours we observed in the sporadic case 
are rather different from those one previously studied [7,8]: 

a) one has a complete synchronization, i.e. l c = N; 

b) this complete synchronization does not change by the addition of a small noise, the 
effect of the noise is a small decreasing of the fraction, r, of the time in which l c = N; 

c) there is a strong correlation between the spatial and temporal behaviour, in 
particular one observes that the effective Lyapunov exponent XAt(t) is large in the 
intervals of time in which l c (t) « AT, while XAt(t) — if l c (t) = N. 
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Figure captions 



Figure 1. Lyapunov exponent A (crosses) and probability of complete 

synchronization r (squares) vs 72. N = 200 and 71 = 0.7. The integration steps 
are 3 • 10 8 for 72 = 0.3, 4.7 ■ 10 s for 72 = 0.2 and a few millions of steps in the other 
cases. 

Figure 2. Time evolution of the mean value defining the Lyapunov exponent A, for 
72 = 0.2 (lower curve) and 72 = 0.3. 

Figure 3. Effective Lyapunov exponent XAt and coherence length l c vs t. N = 200, 
7l = 0.7, 72 = 0.3, At = 3 • 10 4 . 
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